Can the viscosity in astrophysical black hole accretion disks be 
close to its string theory bound? 
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String theory and gauge/gravity duality suggest the lower bound of shear viscosity 
(rj) to entropy density (s) for any matter to be ~ fih/4nkB, when h and ks are 
reduced Planck and Boltzmann constants respectively and \i < 1. Motivated by this, 
we explore r//s in black hole accretion flows, in order to understand if such exotic 
flows could be a natural site for the lowest rj/s. Accretion flow plays an important role 
in black hole physics in identifying the existence of the underlying black hole. This 
is a rotating shear flow with insignificant molecular viscosity, which could however 
have a significant turbulent viscosity, generating transport, heat and hence entropy 
in the flow. However, in presence of strong magnetic field, magnetic stresses can 
help in transporting matter independent of viscosity, via celebrated Blandford-Payne 
mechanism. In such cases, energy and then entropy produces via Ohmic dissipation. 
In addition, certain optically thin, hot, accretion flows, of temperature >10 9 K, may 
be favourable for nuclear burning which could generate/absorb huge energy, much 
higher than that in a star. We find that rj/s in accretion flows appears to be close 
to the lower bound suggested by theory, if they are embedded by strong magnetic 
field or producing nuclear energy, when the source of energy is not viscous effects. 
A lower bound on r//s also leads to an upper bound on the Reynolds number of the 
flow. 
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I. INTRODUCTION 

In order to explain the Quark Gluon Plasma (QGP) observables at Relativistic Heavy Ion 
Collider (RHIC) having temperature (T) >10 12 K, it has become apparent that a non-zero 
shear viscosity, 77, is needed. One way of characterizing 77, which is a dimensionful number, 
is to take its ratio with the entropy density, s. It turns out that RHIC plasma has 

2-0.1*. (i) 

s k B 

Although T) itself is large (in CGS units ~ 10 12 gm/cm/s), its ratio with s is small. Note that 
this ratio for a known fluid, say water, is about two orders of magnitude greater than the 
above value. Hydrodynamic simulations at RHIC and Large Hadron Collider (LHC) (which 
provide the environments created around a few microseconds after the big bang) as well 
suggest a small value for this ratio in the hot QGP state. Such a small value points towards 
strongly interacting matter. In fact, string theory arguments and the famous gauge/gravity 
duality in the form of the AdS/CFT correspondence suggest a lower bound |l| for 77/s for 
any matter given by 

77 u h h 

' > f7-~0.08/ J — , 2 



n nrf n 

1 [2] but non-zero [3|, |4j (see, e.g., [5]). 



where fi < 1 (2| but non-zero [3|, |4j (see, e.g., [5|). For the present purpose, we assume that 
there is a bound and that \i ~ 1. 

Is there any natural site revealing an rj/s close to above lower bound? Note that evenif 
\l « 1, still the same question remains. Such a value may be possible to arise naturally, 
if the temperature and/or density of the systems are/is same as that in RHIC/LHC. In the 
search of such an exotic system, our story starts which could serve as a natural verification 
of gauge/gravity duality in the form of the AdS/CFT correspondence. 

Observed data suggest that the accretion flows around certain black holes must be opti- 
cally thin, geometrically thick and very hot. Such flows often exhibit outflows/jets in their 
hard spectral states. The ion temperature (Tj) in such flows/disks is as high as 7 x 10 11 K. 
The supermassive black hole system at the centre of our galaxy, Sgr A* [6j, (most of the 
temporal classes of) the X-ray binary Cyg X-l [7J are some of the examples of this kind of 
hot accretion disks, which exhibit radiatively inefficient flows. They are particularly differ- 
ent from their radiatively efficient counter part, namely the Keplerian accretion disk, which 



is geometrically thin, optically thick, cooler, and of temperature T ~ 10 K 



a. 



The flow in 



3 

certain temporal classes (and soft spectral states) of micro-quasar GRS 1915+105 [9] is an 
example of Keplerian disk. While the temperature of the latter cases could be same order 
in magnitude as that in the center of a star, due to much low density they practically do 
not exhibit any thermonuclear burning. Hence any source of energy therein must be due 
to magnetic and viscous effects. On the other hand, the former cases may be favourable 
or thermonuclear reactions due to their very high temperature, evenif their density is low 
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. Hence, the energy released/absorbed due to nuclear reactions in an optically thin, 
lot, accretion flow could be comparable to or even dominating over its viscous counter part 
lol . [ill . Q. Being very hot, such flows are also highly ionized and hence expected to be 
strongly magnetized (see e.g. 12]) rendering magnetic dissipation. In fact, the celebrated 
Blandford-Payne mechanism 13]) is based on such a magnetized accretion disk in the Kep- 
lerian regime, which argues for angular momentum transport due to outflows/jets through 
the outgoing magnetic field lines, in absence of shear (turbulent) viscosity. Now an obvious 
question arises: what is the value of nls in optically thin accretion flows (e.g. Advection 
Dominated Accretion Flow (ADAF) 15]. General Advective Accretion Flow (GAAF) fl6|) 
when the temperature is close to that of the QGP matter in RHIC? Primarily guess is that 
it is not too much different from its lower bound as in equation ([1]). Is this that simple to an- 
ticipate? Note that accretion flows with very small molecular viscosity (and then very large 
Reynolds number) exhibit turbulence and then turbulent viscosity. If the turbulent viscos- 
ity dominates, then the energy dissipated due to viscosity in disks is comparable/dominant 
to/over the magnetic/nuclear energy, as will be discussed in detail below. Note also that at 
low densities (zero chemical potential), lattice QCD calculations |l7H19] suggest a crossover 
from the hadronic state to the QGP state at around 1.52 x 10 12 K 2o|. At such high tem- 
peratures, there is expected to be copious pion production in a strongly interacting system. 
This QGP state is thought to have existed a few microseconds after the big bang and is 
being extensively studied at RHIC and LHC. Of course the density in early universe is large 
as well, unlike that in accretion flows, a large temperature is well enough to reveal such a 
phase which has been mimicked in RHIC/LHC. However, at a very large density (several 
orders of magnitude larger than that in accretion flows) the lattice QCD results are not valid 
(Zoltan Fodor, private communication) . 



Recently, Sinha & Mukhopadhyay 



21] initiated to look at the above issue and argued that 



the Shakura-Sunyaev turbulent viscosity parameter a [8] should not be constant throughout 
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the flow, and be decreasing with the increase of temperature and/or density of the flow. 
Based on the general relativistic model of a viscous ADAF, they showed that rj/s for an 
astrophysical black hole, at finite a, is always several orders of magnitude higher than its 
value in a QGP fluid. Then they argued for the flow to become very weekly viscous (not 
exhibiting turbulence) close to a black hole, rendering a smaller rj/s. The results appear to 
be independent of the choice of equation of state (EoS); whether of ideal gas or of QCD which 
is appropriate for flows with T>10 12 K and low density (close to zero chemical potential). 
However, they could not clarify the natural circumstances when a is small. Note that small 
a would reveal negligible dissipation of energy in the flow. Hence, question arises what is 
the source of energy and entropy in the accretion flows, which is also very important in 
order to address 77/s? In addition, how is the accretion possible in a very small-a flow? 
Moreover, as will be explained in the next section, above work did not model the entropy of 
the flow adequately, while generally interpreted the results correctly. Instead, that was the 
first step forward in the direction of evaluating 77/s for any astrophysical flow, eventhough 
incomplete, which is a potential natural site for small rj/s. In the present work, we have 
removed the above uncertainties lying in the previous work. Then we plan to establish the 
underlying physics giving rise to the astrophysical values of rj/s. More precisely, we plan 
to address: (1) regime of accretion flows giving rise to rj/s close to its theoretical lower 
bound, (2) observational implications of such regimes, (3) black hole sources presumably 
revealing lower bound of rj/s, (4) constraining physical parameters (e.g. Reynolds number) 
of accretion flows based on rj/s bound. 

In the next section, we discuss the salient features of hot advective accretion flows and 
the importance of strong magnetic field and nuclear energy produced therein in order to 
obtain small rj/s. Subsequently in §3 we analyse the solutions of the equations and address 
the flow parameters making rj/s to be close to the lower limit. Finally §4 summarizes our 
findings with discussion and implication. 

II. ADVECTIVE ACCRETION FLOWS AROUND BLACK HOLES 



In order to understand the underlying physics, let us consider the equations describing 
optically thin, viscous, magnetized, advective accretion disk {24]. We describe the model in 
the pseudo-Newtonian framework with the Paczyhski-Wiita potential 25|, |26[. This is par- 
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ticularly because, for the present purpose, the pure general relativistic results for a rotating 
black hole qualitatively do not reveal any new physics. Hence, the vertically averaged hy- 
dromagnetic equations of energy-momentum balance along with the equations of continuity, 
induction and divergence of magnetic field in the limit of very large conductivity are 

M = 4nxm, 

^dd IdP _ A 2 1 _ ( , 

dx p dx x 3 2(x — l) 2 

S dX 1 d (V W \ ,hx ( dB^ dB^ B X B. 
v— - [x W X(j) ) + — B x — h B z — 1 

dx xh dx L \ dx dz x 

dx T3 — 1 \ dx p dx J 

± (xSB z ) = 0, (3) 

assuming that the variables do not vary significantly in the vertical direction such that 
d/dz — > 1/z, which is indeed true for the disk flows. Here M is the conserved mass accretion 
rate, E and p are the vertically integrated density and density of the flow respectively, 
d is the radial velocity, P the total pressure, A the specific angular momentum, W x <f> the 
vertically integrated shearing stress, h the half-thickness, s the entropy per unit volume, T 
the temperature of the flow, Q + and Q~ are the vertically integrated net energy released and 
absorbed rates in/from the flow respectively, r 1; T 3 indicate the polytropic indices depending 
on the gas and radiation content in the flow (see, e.g., [16| for exact expressions) and B x , B^ 
and B z are the components of magnetic field. Note that, the independent variable x is the 
radial coordinate of the flow expressed in the units of 2GM/c 2 , where G is the gravitation 
constant, M the mass of the black hole and c the speed of light. Accordingly, all the 
above variab 
literature 

Let us now pay a special attention to the energy equation of the equation set ([3]). We 
can rewrite it as 

ds 

ax = ^ vis ^ ^ mag ^ ^ Nex ~ ^ md ~ ^ Nend ' (■ ' 

where g vis , g mag and gNex respectively define the energies released per unit volume per unit 



es are expressed in dimensionless units. For any other details, see the existing 
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time due to viscous dissipation, magnetic dissipation and thermonuclear reactions and g ra d 
and gNend respectively indicate the energy radiated away per unit volume per unit time by 
various cooling processes like Bremsstrahlung, synchrotron and inverse-Comptonization of 
soft photons, and the energy absorbed per unit volume per unit time due to thermonuclear 
reactions. 

It was shown earlier that the nuclear energy released/absorbed in a hot accretion disk 



need not be negligible 
over-abundant metals 



10 



10 



29 



4. 



and the thermonuclear reactions therein could produce 



3l| which are observed in galaxies. However, the value of a 



viscosity determines whether the energy due to nuclear reactions is comparable to that due 



to the viscous effects 
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and Ohmic dissipation or not. From equation fljj) we can write 
down the entropy per unit volume, namely the entropy density, of the flow as 

Qvis + Qmag + ^Nex — Qrad — (/Nend 



■ dx. 



(5) 



Now the turbulent kinematic viscosity v = $ t l t /3, when *& t and l t are respectively the 
turbulent Eddy velocity and the size of Eddy, which are respectively to be scaled linearly 
with sound speed (c s ) of the disk and h and hence 

a$c s a t h 



ac s h, 



where a#, on and hence a are constants. Therefore, the shear viscosity is given by 



r] = apc s h. 



(6) 



(7) 



In an optically thin flow, often g rad is approximately chosen to be proportional to g vis + 
<Zma g , in order to understand the hydrodynamic properties of the flow without much loss of 
generality, such that q ra & = (1 — /)(<? v is + <7mag) (see, e.g., Q, Q); when / is a constant. 
For an ADAF / = 1 strictly (and M « M Edd , when Msdd is the Eddington accretion 
rate) and for a GAAF < / < 1 (and M is larger than that in an ADAF |l6j). Note that 
f = for an optically thick Keplerian accretion disk. Hence, following previous work (e.g. 



IS 



24 



271]). we can write 



<?vis + gmag - grad = / (<?vis + ^mag) = / ( a(/n+l^ + ^^P^ ( ~ J + ^ J > ( 8 ) 

where I n = (2 n n\) 2 / (2n + 1)!. Therefore, if g mag , quejf, q^end/ f « <?vis, which we name 
as viscous regime, rj/s does not explicitly depend on a (and 77) [37|, depends only on the 
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hydro dynamic quantities whose numerical values, in a fixed range of radii, do not vary 
significantly in the optically thin regime with the change of a. In the viscous regime, rj/s 
simply scales with M. Note that viscosity depends on the size of the system, which in turn 
is determined by the mass of the black hole. Hence, only for a primordial black hole of 
M<10 16 gm, rj/s is close to its theoretical lower bound, as was suggested earlier [2l|. It was 
then also suggested that a has to be very small in an ADAF to obtain rj/s ~ h/AirkB for 
an astrophysical black hole. However, for a small a, viscous dissipation (and then s) in an 
ADAF (or generally in a GAAF) should be small as well, as is apparent from equations (j4]) 
and (JED- Hence, the ratio of rj to s should not be invariably small, as described above, in 
the viscous regime. Therefore, while the earlier work 2l| is very interesting and the first 
attempt to evaluate 77/5 in an astrophysical flow, the flow (and hence any conclusion) should 
be amended with extra physics. In other words, the conclusion is incomplete. 

Indeed r]/s could be small at a small a when s does not decrease in that extent 
with the decrease of a. However, this is possible only when g mag >> fq v[s or/and 
(/Nex or/and gNcnd >> /<?vis, which we name magnetic and nuclear regimes respectively. The 
magnetic regime essentially corresponds to the Blandford-Payne mechanism in the hard 
spectral state, when accretion takes place via magnetic stress, independent of the presence 
of shear/turbulent viscosity [3]. This plausibly corresponds to the case of accretion flows 
with outflows/jets, e.g. Sgr A*, hard states of GRS 1915+105. Hence, even at a negligibly 
small a accretion is possible with a negligible rj/s. On the other hand, the nuclear reactions 
and correspondin g en ergi es depend on the density and temperature of the system. As shown 
32M34J] that even in the inviscid limit, temperature and density of the op- 



earlier 
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tically thin flows remain as high as those of a viscous accretion flow for a fixed M, which 



successfully explain many observation [35|, [36( . Indeed, it was shown lOj, [ljj that significant 
nuclear energy (~ 10 15 ergs/cc/sec) throughout the disk may be released/absorbed in an 
inviscid flow which could lead to disk instability. Hence, from equation ([5l) it is naturally 
possible that s is large even at a small a. Therefore, rj/s ~ h/AirkB also plausibly in the 



nuclear regime of the accretion flow, when 

m 



the case in the hot accretion flows 
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JffNe 

1 



<?Nend| ~ |<?Ne X | or |g Ne nd| which could be 
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III. SOLUTIONS AND Tj / S 

Let us understand the above physics quantitatively. Figures [TK,b,c show hydrody- 
namic/hydromagnetic properties of a disk around a 1OM nonrotating black hole in the 
magnetic regime, accreting mass at the rate ~ ME dd . in this case of weakly viscous hot flow 
with q~5x 10~ 17 (but not zero exactly), rj/s is shown to be close to the theoretical lower 
bound h/4irkB- For comparison, rj/s is also shown for a viscous regime, which is quite large 
for a = 0.05 and / = 0.1. We also show the variations of Mach number and magnetic field in 
the disk for the cases mentioned above for reference. Interestingly, increase of magnetic field 
in the magnetic regime increases rate of transport which further renders advancing the Ke- 
plerian disk towards the black hole, making the sub-Keplerian regime, which is shown here, 
smaller. Figures HH,e,f depict the properties for M ~ 10 _2 ME dd i n the magnetic regime. 
A very small (but nonzero) a (~ 2 x 10 -19 ) reveals r]/s to exhibit the theoretical lower 
limit. The impact of change of the magnetic field strength to the disk is same as that in the 
high accreting case discussed above. Interestingly, independent of the accretion rate, mag- 
netic field >kG is able to transport matter efficiently in the sub-Keplerian flow in absence 



of viscosity, as was already argued by Blandford & Payne [13( in the Keplerian, self-similar 
framework. A smaller field strength would also transport matter as well, except with the 
increase of size of the sub-Keplerian regime. 

Figure [2] is devoted for hydrodynamic/hydromagnetic properties around a 1O 7 M non- 
rotating black hole in the magnetic regime. While Figs. EK,b,c correspond to M ~ M^dd, 
Figs. [2]i,e,f correspond to M ~ 10 _4 MEdd- In either of the cases, very weakly viscous 
(a ~ 5 x 10~ 23 for M ~ M Edd and a ~ 10~ 25 for M ~ 10" 4 M Edd ) hot flows reveal small 
rj/s, close to its theoretical lower bound h/4nkB- For comparison, rj/s is also shown for a 
viscous, high accreting regime which is quite large for a = 0.05 and / = 0.5. Interestingly, 
magnetic field >mG is able to transport matter efficiently in the sub-Keplerian flow around 
a supermassive black hole in absence of viscosity, as was already argued by Blandford & 



Payne 13( in the Keplerian, self-similar framework. As is the case for a stellar mass black 
hole, a smaller field strength would also transport matter, except the fact that the Keplerian 
to sub-Keplerian transition would arise at a larger radius. 

Let us also look into the nuclear regime. Figure |3] shows hydrodynamic properties of a 
disk around a 1OM nonrotating black with M = 0.5MEdd in the nuclear regime. In the cases 
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FIG. 1: Variations of (a) Mach number, (b) magnetic field in G, (c) rj/s in CGS unit when 
the horizontal line (long-dashed) refers to h/inks, as functions of flow radius for M = WMq 
at M = MEdd- Solid line corresponds to the flow in viscous regime and dotted and dashed lines 
correspond to the flows in magnetic regime, (d), (e), (f) describe same as in (a), (b), (c) respectively, 
except for M = O.OlMEdd- See Table 1 for other details. 

of weakly viscous hot flow, rj/s is shown to be inversely proportional to the nuclear energy 
released per unit volume. When the ratio of nuclear energy per unit volume to viscous energy 
per unit volume is chosen to be Qn/QvIs ~ 10 16 (which corresponds to Q ma , x ~ 3.4 x 10 17 
erg/sec, depicted by the solid line), rj/s appears to be close to the theoretical lower bound 
h/inks- Note that the production of nuclear energy also depends on the residence time of 
matter in the disk which is determined by a, polytropic constant 7, initial mass fraction of 
the inflowing isotopes etc. For comparison, rj/s is also shown for a viscous regime, which 
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FIG. 2: Variations of (a) Mach number, (b) magnetic field in G, (c) r//s in CGS unit when the 
horizontal line (long-dashed) refers to h/AnkB, as functions of flow radius for M = 1O 7 M at 
M = MEdd- Solid line corresponds to the flow in viscous regime and dotted line corresponds to 
the flow in magnetic regime, (d), (e), (f) describe same as in (a), (b), (c) respectively, except for 
M = 10" 4 M E dd- See Table 2 for other details. 

is quite large for a = 0.05 and / = 0.1. We also show the variations of Mach number 
and temperature in the disk for the cases mentioned above for reference. Interestingly, the 
latter profiles overlap each other in the nuclear regime, independent of Q^/Qvis, and are 
very similar to that in the viscous regime. This also verifies that the Mach number and 
temperature do not depend much on a, in the optically thin flows. 

In Fig. HI we show a couple of sets of solutions for highly sub-Eddington hotter flows and 
a solution set for a low angular momentum, high mass accreting, hot flow, in the nuclear 
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regimes. Interestingly, a very low angular momentum (A ~ 1), gas dominated, hot accretion 
flow with M = 10 _4 MEdd (solid line) violates the theoretical lower limit of rj/s. This restricts 
A of the flow. If A increases from 1 to 2, the flow becomes relatively cooler due to decrease of 
radial velocity which restricts nuclear energy, rendering slightly larger 77/s > h/AnkB (dotted 
line). A flow with similar A (~ 1.6) and a, but M = 0.5MEdd and hence revealing a radiation 
dominated accretion (dashed line) flow, exhibits however a similar rj/s > h/4irkB- Note that 
while the nuclear energy among the three cases considered here is highest in the last flow, 
due to presence of radiation the cooling efficiency therein is maximum, which renders its 77/s 
above the lower limit throughout. 

Table 1 

Parameters for the results shown in Figure [T] when M = 1OM 



M 


7 


A 


Regime 




\B\ in CGS 


/ 


Figure, line 


in Msdd units 




in inner disk 






in inner disk 






1 


1.335 


3.2 


Magnetic, a = 5 x 10~ 


-17 


1.85 x 10 3 G 


0.5 


Figs. [TK,b,c, dotted 


1 


1.335 


3.2 


Magnetic, a = 5 x 10" 


-17 


1.85 x 10 5 G 


0.5 


Figs. [T^,b,c, dashed 


1 


1.335 


3.2 


Viscous, a = 0.05 






0.1 


Figs. [Tk,c, solid 


0.01 


1.5 


3.25 


Magnetic, a = 2 x 10~ 


-19 


220G 


0.5 


Figs. HH,e,f, dotted 


0.01 


1.5 


3.25 


Magnetic, a = 2 x 10~ 


-19 


2.2 x 10 4 G 


0.5 


Figs. [T]i,e,f, dashed 



Table 2 

Parameters for the results shown in Figure [2] when M = 1O 7 M 



M 


7 


A 


Reg 


ime 


\B\ in CGS 


/ 


Figure, line 


in MEdd units 




in inner disk 






in inner disk 
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1.335 


3.2 


Magnetic, a 


= 5 x 10~ 23 


1.85G 


0.5 


Figs. [2k,b,c, dotted 


1 


1.335 


3.2 


Viscous, 


a = 0.05 




0.5 


Figs. |2^,c, solid 


10" 4 


1.53 


2.6 


Magnetic, 


a = 10~ 25 


0.0185G 


0.5 


Figs. EH,e,f, dotted 
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FIG. 3: Variations of (a) Mach number, (b) temperature in K, (c) energy density in CGS unit, (d) 
rf/s in CGS unit, as functions of flow radius for M = 10M Q at high mass accretion rates. Solid 
and dotted lines correspond to the flows of very high and high nuclear regimes respectively and 
dashed line corresponds to the flow of viscous regime. See Table 3 for other details. 
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Table 3 



Parameters for the results shown in Figures El and H] when M = 10M© 



M in MEdd units 


7 


A 


Regime 


Qmax in CGS 
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Fi| 


jure, line 


0.5 


1.34 


3.2 


Nuclear 




3.4 x 10 17 




Fig 


;. El solid 
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1.34 


3.2 


Nuclear 




3.4 x 10 11 




Fig. 


El dotted 
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1.335 


3.2 (inner) 


Viscous, a = 


0.05 


1.8 x 10 19 


0.1 


Fig. 


El dashed 


0.0001 


1.53 
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Nuclear 




1.3 x 10 14 




Fig 


;. HI solid 


0.0001 


1.53 


2 


Nuclear 




3 x 10 12 




Fig. 


H dotted 


0.5 


1.34 


1.6 


Nuclear 




8.4 x 10 19 




Fig. 


IH dashed 



IV. SUMMARY AND DISCUSSION 



We have shown the astrophysical existence of small T]/s, close to its theoretical lower 
bound. As our physical inputs, governing the values of rj/s, are in accordance with observed 
data, this could be argued as a natural evidence for r\/s close to its lower bound. This is 
possible in the cases of hot, optically thin, accretion flows. More precisely, in order to have 
a small 77/s, flow should be magnetically dominated and/or producing huge energy /entropy 
other than that due to viscosity. 

The celebrated Blandford-Payne mechanism [l3| was already shown to produce accretion 
via magnetic stress (which helps in producing outflows/jets and removing angular momen- 
tum) in absence of viscosity in the framework of Keplerian, self-similar flow. We have 
employed this idea, but relaxing self-similar assumption, in order to obtain 77/s close to 
h/AnkB in our accretion flows. We have found that in the presence of magnetic field of 
the order of mG to kG, the disk flows exhibit this value of rj/s for supermassive to stellar 
mass black holes. Observed hard X-rays in Cyg X-l, some temporal classes (e.g. x) °f 
GRS 1915+105, Sgr A* are observed to be associated with jets which argues for the sources 
to exhibit strong magnetic fields. Hence, those sources could be considered to be natural 
sites for 77/s ~ h/AirkB- 

Additionally, the flow exhibiting nuclear energy could be a natural site for small 77/s. 
Nuclear energy depends on the density and temperature of the flow. The temperature in 
an optically thin, sub-Keplerian accretion flow, producing observed hard X-rays, could be 



14 




FIG. 4: Variations of (a) Mach number, (b) temperature in K, (c) energy density in CGS unit, (d) 
rj/s in CGS unit, as functions of flow radius for M = WMq. Solid and dotted lines correspond 
to the flows of very high and high nuclear regimes respectively at low mass accretion rates and 
dashed line corresponds to the flow of a nuclear regime at a high mass accretion rate. See Table 3 
for other details. 

>10 9 K. Therefore, eventhough the density therein is small, very high temperature would 
suffice for nuclear burning. We have recalled the previous work discussing pos sible nuclear 



burning in accretion flows and its observational evidences/consequences jiol . [ill H-Sl]. 
Based on them, we have modeled the hot, optically thin accretion flows for the present 
purpose and shown how the nuclear energy in the flow determines the value of r]/s. Thus 
we argue that the flows producing hard X-rays in Cyg X-l, some temporal classes (e.g. 
x) of GRS 1915+105 could be hot enough to have rj/s ~ h/4iik B , atleast close to the 
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black hole for certain nuclear regimes, when flows are necessarily sub-Keplerian due to the 
dominance of gravitational power of the black hole. It will be now interesting to explore rj/s 
for other astrophysical flows, having, e.g., strong magnetic field or/and high temperature 
and/or density. 

Finally, we end by enlightening maximum Reynolds number (R e ) to the flows, restricted 
by the lower bound of rj/s. R e , which is the ratio of typical (characteristic) velocity times 
length scale to the kinematic viscosity of the flow, can be recast as R e = M ac (x/h)/a. Now 
from Figs. [1] and [3j M ac close to the black hole is ~ 5 — 10. However, the magnetic and 
nuclear regimes correspond to very low but nonzero viscosity revealing a>10 -16 for rj/s lower 
bound to be satisfied. As typically h/x ~ 0.1 — 0.5 for an optically thin, geometrically thick 
flow, for stellar mass black holes i? e <10 17 . If we look at Fig. [2] showing M ac to be of the 
similar order in magnitude around supermassive black holes, however for smaller values of a 
(>10 -22 ) in larger accretion disks (as size scales with the mass of the black hole), we obtain 
i? e <10 23 for those flows. 
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